Block-block entanglement and quantum phase transitions in one-dimensional 
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In this paper, we study block-block entanglement in the ground state of one-dimensional ex- 
tended Hubbard model. Our results show that the phase diagram derived from the block-block 
entanglement manifests richer structure than that of the local (single site) entanglement because it 
comprises nonlocal correlation. Besides phases characterized by the charge-density-wave, the spin- 
density-wave, and phase-separation, which can be sketched out by the local entanglement, singlet 
superconductivity phase could be identified on the contour map of the block-block entanglement. 
Scaling analysis shows that logj (l) behavior of the block-block entanglement may exist in both non- 
critical and the critical regions, while some local extremum are induced by the finite-size effect. We 
also study the block-block entanglement defined in the momentum space and discuss its relation to 
the phase transition from singlet superconducting state to the charge-density-wave state. 
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I. INTRODUCTION 

Entanglement, as one of the most intriguing feature 
of quantum mechanics^ and a crucial resource in quan- 
tum information theor}i2i^, has been a subject of much 
study in recent years. For the purpose of practical appli- 
cation, to realize quantum information processing, such 
as quantum state transfer*- , based on the condensed 
matter physics is of great importance. Such potential 
prospect in the application now has motivated many the- 
oretical investigations on the entanglement in spin and 
fermionic systemsSiSi^iiii. On the other hand, quan- 
tum phase transitions (QPTs), which are characterized 
by change in the properties of the ground state of a 
many-body system, are generally driven by quantum 
fluctuation at zero temperature. One therefore expects 
that the entanglement, as a term of pure quantum cor- 
relation, should closely relate to QPTs. Indeed, some 
observed relationshipo"ii2ii3ii4,i5,i6,i7,i8,i9,20,2i between 

the entanglement and QFTs^S now have not only deep- 
ened our understanding on the QPTs, but also strength- 
ened the connection between the quantum information 
theory and condensed matter physics. For example, Os- 
terloh et al.j^ reported that the pairwise entanglement of 
two nearest neighbors shows scaling behavior in the vicin- 
ity of quantum phase transition point of the transverse- 
field Ising model. For other models, such as the XXZ 
modeii^, spin model with mutual exchangesi^, etc., the 
pairwise entanglement also manifests various interesting 
properties, such as being maximum at transition point, 
exhibits singularity, and shows scaling behavior. 

Besides the pairwise entanglement between two sites, 
for a non-degenerate ground state, the block-block 
entanglemen^iSiiS also provides a good quantity to de- 
scribe the pure quantum correlation at zero temperature. 
The scaling property of the block-block entanglement 
establishes an interesting connection between concepts 
of quantum information and quantum field theorySS,. 



Moreover, unlike the concurrence^^, a measurement of 
the entanglement of formation, which is computable 
only for spin- 1/2 system, the block-block entanglement 
can be generalized to high-spin and fermionic systems. 
There are some works which studied the entanglement 
in fermionic systems2ii2Si2&. However, the investigations 
on the entanglement in relation to QPTs are still on an 
early stago^7i2?i2?i?9i?i. 

In the extended Hubbard model, it has been shown 
that the global phase diagram can be sketched out 
by the contour map of a quantity called the local 
entanglementSSiSl. However, the local entanglement fails 
to identify phases which are related to the off-diagonal- 
long-range order, such as the superconducting phase. In 
this paper, we extend our previous investigation on the 
local entanglement in the extended Hubbard model to 
the block-block entanglement. We show that the block- 
block entanglement with block size larger than one can 
provide more useful information than the local entangle- 
ment since it comprises the nonlocal correlation in its 
expression. The paper is organized as follows. In Sec. ^ 
we first briefly review some basic knowledge of the ex- 
tended Hubbard model and then introduce the definition 
of the block-block entanglement. In Sec. IIIII we address 
the problem of the global phase diagram on the U — V 
plane by introducing the contour map of the block-block 
entanglement. In Sec. IIVI we study the scaling behavior 
of the block-block entanglement by changing the system 
size. One of the interesting results obtained in this work 
is that the log2(^) behavior {I is the block size) of the 
block-block entanglement also shows in noncritical re- 
gion. In Sec. EI we study the finite-size effort and clarify 
some unaccustomed behaviors of the block-block entan- 
glement. A simple expression of the local entanglement 
is obtained. In Sec. IVII we try other measure of the 
block-block entanglement in order to explore the phase 
boundary which may not be identified by the extremum 
of the local entanglement. Finally, a brief summary is 
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given in Sec. IVIII 

II. THE MODEL AND FORMULISM 

The one-dimensional extended Hubbard model is de- 
fined by the Hamiltonian 

o-,i:i5 j 3 

where a =t,i;i = 5 = ±1, c j ^ and Cj^a are 

creation and annihilation operators at site j, U and V 
define the on-site and nearest-neighbor Coulomb inter- 
actions, respectively. On each site the local states have 
four possible configurations, denoted by 

1,2,3,4)) = |0), IT), U), lU). (2) 

The Hilbert space associated with L-site system, known 
as the Fock space Hf{L), is spanned by 4^ basis vec- 
tors \ji,...,jL) = rij'^i l<^j(^i)>- Any state in such 
a system can be expressed as a superposition of these 
basis. If we choose periodic boundary condition for 
L = An + 2 and anti-periodic boundary condition for 
L = An, where n is an arbitrary integer, the ground state 
is non-degenerate^. Therefore, considering the reduced 
density matrix of a block of I successive site of the ground 
state pi = trr|^)(4'|, the von Neumann entropy £'i,(Z), i.e. 

E^il) = -tT[pi\og^{pi)] (3) 

measures the entanglement between the / sites and the 
rest L — I sites of the system, as shown in the following 
figure for 10-site system with I — 4. 

A simple case of the block-block entanglement is when 
only one site is taken into account. Then the reduced 
density matrix can be written into a simple form.^S*2X, 

Pi^z\ o)(o I +U+ |T)(TI +u- +w |Ti)(Til, (4) 

in which 

1^ = ("-T"-i> = Tr(n|n|pi), 

— (rif) — w, = (rij) — w, 
z — \ ~ ^ — w ~ 1 — (n|) — {ui) + w. (5) 

Here (n^) and {n^) are electron densities with spin up and 
spin down respectively. The corresponding von Neumann 
entropy then takes the form 

El, = —Z log2 Z — log2 — U" log2 U~ — W log2 W . 

If / > 1, the reduced density matrix can not be written 
into a simple expression. However, since the Hamiltonian 



(P) has U(l)xSU(2) symmetry: Cj^a e'^Cj,^; Cj,o- ^ 
UaSCj^Si which manifests the charge conservation and in- 
variance under spin rotation UaS- Therefore, any eigen- 
state of the Hamiltonian iQJ can be both the eigenstate of 
z component of total spins and of particle number. This 
fact leads to that, for arbitrary block size I, there is no 
coherent superposition of local state with different value 
of 5^(0 = Ej-K.T -^j.i) and N(l) = Ej K',? + '^jm) in 
the reduced density matrix. That is, the reduced density 
matrix must have the diagonal form classified by both 
S^{1) and N{1). In this paper, we apply exact diagonal- 
ization technique to get the ground state and then exactly 
diagonalize the reduced density matrix to calculate the 
block-block entanglement. 



III. ENTANGLEMENT AND THE PHASE 
DIAGRAM 

The extended Hubbard model is a prototype model 
in condensed matter theory for it exhibits a rich phase 
diagram22i2li2^ where various quantum phase transitions 
occur among symmetry broken states. The correspond- 
ing symmetry broken states typically include the charge- 
density- wave (CDW), the spin-density- wave (SDW), 
phase separation (PS), singlet (SS) and triplet supercon- 
ducting phase(TS), and bond-order wave (B0W)^Si2LS. 
Many efforts have been made to obtain the phase dia- 
gram of the extended Hubbard model at different band 
fillings on the U-V plane^^i2Si^L2&. In our previous work 
on the local entanglement21, it is remarkable to see that 
the skeleton of the phase diagram of the extended Hub- 
bard model can be directly obtained from the contour 
map of the local entanglement. This is by no means triv- 
ial. In the conventional approach to obtain the phase di- 
agram of the extended Hubbard model, one has to study 
behaviors of different order parameters in different re- 
gions, either by comparing the ground-state energy or the 
critical exponent of correlation functions associated with 
broken symmetry. While using a single quantity, , the 
global picture of the system at zero temperature can be 
observed. This is not a coincident, for the non- vanishing 
order parameter means the existence of non-local corre- 
sponding correlation, whose quantum part is just the en- 
tanglement, so the competitions between different orders 
may lead to changes in the entanglement. When a QPT 
occurs, entanglement will also behave distinctively in dif- 
ferent phases. Therefore, this result reflects the underly- 
ing relation between the entanglement and QPTs beyond 
the superposition principle of quantum mechanics. How- 
ever, there are some limitations. For example, the local 
entanglement can not be used to identify superconduct- 
ing phases, due to the fact that the broken symmetry 
is associated with off-diagonal-long-range-order, whereas 
the local reduced density matrix is diagonal. To identify 
such phase is a challenge, and it is one of our motivations 
to include the off-diagonal correlation in the study of the 
entanglement and QPTs. 
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FIG. 1: The block-block entanglement and its contour map 



In Fig. ^ we show a three-dimensional diagram, as 
well as its contour map of the block-block entanglement 
for 8 sites system with block size ^ = 3. Obviously 
the structure of the contour map is richer than that of 
the local entanglement-2i, especially at the region near 
U <0, — l<y<0. The contour lines should be simi- 
lar in the same phase region, because the entanglement 
should behave in a similar way in the same phase. So 
the distinct change of the contour lines suggests the ex- 
istence of a new phase. Compare this contour map with 
the phase diagram of the extended Hubbard model at 
half-filling, this region corresponds to SS phase. For the 
transition from CDW to SDW, the block-block entangle- 
ment, including the local entanglement^'' exhibits maxi- 
mal value. This fact can be understood by noting that, at 
the QPT point from CDW to SDW, the weights for differ- 
ent component in the reduced-density matrix are closer 
due to the higher symmetry in the critical region. Take 
the local entanglement as an example, we have equally 
valued and z at the transition point. From 

this point of view, the extremum behavior of the entan- 
glement at the transition point manifests that the QPTs 
are induced by symmetry broken^^. For the transition 
from SS to CDW, the contour lines have an inflexion, 
which suggests the existence of other types of phase tran- 
sitions. 
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FIG. 2: The block-block entanglement as a function of V for 
various block size. Here L — 10,U — —2. 



IV. ENTANGLEMENT CHANGE WITH BLOCK 

SIZE 

In Fig. 12 we show the block-block entanglement as 
a function of V for various block size at fixed U = —2, 
L = 10. In the cross section of U = —2, there are three 
phases. If y <C — 1, the ground state is phase separated. 
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FIG. 3: The block-block entanglement as a function of V for 
various block size. Here L = 10,U = 0. 



That is almost half of sites are doubly occupied and con- 
gregate together, while another half of sites are empty. 
The corresponding wavefunction is dominated by the fol- 
lowing configuration, e.g. for the 10-sites system, 

PS(a) : Ti Ti Ti Ti Ti 0, 

which is in fact an eigenstate of the Hamiltonian with 
t = 0. Obviously, this state is separable and has no en- 
tanglement. However, in the presence of the hoping term, 
the ground state becomes a superposition of all possible 
configurations transformed from the above one through 
an arbitrary translation operation. This process intro- 
duces entanglement into the ground state. If we take the 
reduced-density matrix with block size I, the number of 
block matrix corresponding to different particle numbers, 
N{1) = 0, 2, • • • , 2/ is / + 1. The dimension of the block 
matrix with a given nonzero N{1) — 2n, n < is l — n + 1, 
then the number of nonzero eigenvalues of the reduced 
density matrix pi is {P + I + 2)/2. This result implies 
a log(Z) behavior of the block-block entanglement in the 
phase-separated region. That is, the block-block entan- 
glement will not tend to a constant with the increasing of 
block size, as shown in Fig. |2| This phenomenon is quite 
different from previous studies of spin chains^^, where 
the log(/) behavior of the block-block entanglement only 
exists at the critical region, while at noncritical region, 
it will tend to a constant as I increases. 

When the absolute value of negative V becomes 
smaller, the effect of the hoping term and the on-site in- 
teraction can not be neglected. Though the later still try 
to maintain the number of doubly occupied sites, the for- 
mer really want to diffuse the cumulated electron pairs. 
Then the configuration PS(a) will not be the dominant 
one, and the weights of the following configurations 

PS(b) : Ti Ti Ti Ti T i 



PS(c) : Ti Ti Ti Ti Ti 



will increase. Obviously the configuration PS(b) will in- 
troduce block matrices with odd number of particle into 
the reduced-density matrix, and thus increase the entan- 
glement. So the ordered phase is destroyed gradually as 

V increases, and finally the system has maximum entan- 
glement and transits to the SS phase. In the SS phase, 
the ground state is characterized by an off-diagonal-long- 
range order. U V ^ 0, the off-diagonal correlation func- 
tions in the reduced-density matrix will cause a decrease 
in entanglement which is consistent with the figure. If 

V > 0, the hoping process will be suppressed, and the 
off-diagonal-long-range order will be destroyed. This fact 
leads to further decrease of the block-block entanglement. 
Moreover, with the increasing of V, the weight of CDW 
will increase almost to 1, so that the reduced density 
matrix is composed of 1/2 Ti and each, and the corre- 
sponding entanglement can be calculated exactly: 



E^l) 



(6) 



1=1,2 



This result is consistent with the previous work on the 
spin modeliS, which says that the block-block entangle- 
ment will tend to a constant in the noncritical region as 
the block size increases. 

In Fig. 13 we show the block-block entanglement as a 
function of V for various block size at fixed U = 0,1 = 
10. Similar analysis for Fig. [3 can be applied to Fig. 
|31 The main difference is that in the nearby region of 
y = 0, it seems to be a flat area. As SDW, CDW, SS 
and TS all exist in this region, it is not easy to divide 
them clearly by using the block-block entanglement only. 
We also take more detailed calculation in this region. It 
is surprising that for L = 10, U = 0, the block-block 
entanglement get a minimum value at V = 0.0, and a 
maximum value at = 0.09. It seems that the minimum 
point is a transition point, just like the appearance of 
a maximum in the local entanglement^ at this point. 
However, whether these points witness phase transitions 
needs further investigation on longer chains. 

^From the above investigations, we find that the block- 
block entanglement is an increasing function of I in the 
region I g [0,L/2]. In our study, in order to obtain the 
entropy of reduced density matrix, we need to diagonalize 
density matrix by the standard QR algorithm. Limited 
by the computer power, the largest size we have studied 
is 12, then the block size varies from 1 to 6. It has been 
argued that the behavior of the entanglement in a crit- 
ical spin chain matches the result of the conformal field 
theory, where the geometry entropy is analogous to the 
spin block entropy. Although we can not give the exact 
result of scaling on the EHM, we can still see from Fig. 0| 
that the block-block entanglement grows logarithmically 
with the block size in some critical regions. Moreover, as 
we pointed out in the above, such a logarithmic behav- 
ior of the block-block entanglement can also exist in the 
noncritical region for the extended Hubbard model. 
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FIG. 5: The block-block entanglement change with different 
values of V for fixed block size I = 4 and various system size. 
Here V = -4. 



ENTANGLEMENT CHANGE WITH 
SYSTEM SIZE 



In this section, we study the dependence of the block- 
block entanglement on the system size. In our previous 
work on the local entanglement of the Hubbard model, 
the local entanglement in some regions (CDW and SDW, 
for example) obtained from the Bethe Ansatz method 
(for L — oo and L = 70) and by the exact diagonaliza- 
tion technique (L = 10) agree with each other excellently. 
However, if one wants to study phases such as phase 
separation, entanglement will show strong finite size de- 
pendence. For simplicity, we look at the local entangle- 
ment again. If U is negative and ^ ^ — 1, the ground 
state is just a superposition of the configurations trans- 
formed from PS(a) by the translation operation. The 
components in pi are z = w ~ 1/2, u+ = ~ 0, then 
Ey{l) ~ 1. However, if the on-site interaction becomes 
positive and larger, the ground-state will prefer the con- 



figuration 

PS(d) :0 



i Ti Ti n Ti T 0, 



to the PS(a). For the state dominated by the PS(d) 
configuration, the local entanglement has a strong finite 
size dependence; 



El 



log2 L 



1 



l0g2 



(7) 



Though it will tend to 1 in the thermodynamic limit, it 
develops a jump in the 3D figure of the local entangle- 
ment for a small systemSi. Such a jump in a finite chain 
signals a crossover from one PS configurations to another. 

We now address the problem for the case of ^ > 1. We 
show the block-block entanglement change with the value 
of V for fixed I = 4 and vary system size in Fig. [S] From 
the figure, we first notice a jump around U — 5, which 
is the critical point for the transition from PS to SDW. 
We also see that there are seemly jumps of the block- 
block entanglement around U = 0.6. However, it is not a 
critical point. It is just a result of finite-size effect, as we 
found in the local entanglement. Because in both U < 0.6 
and 0.6 < C/ < 4 regions, the ground state is dominated 
by phase separation configurations. The difference is that 
in the region U < 0.6, the dominant configuration is 
PS (a). As U increases, one doubly occupied site tends to 
be singly occupied. So the configuration PS(d) becomes 
the dominating configuration. Therefore the jump here 
does not correspond to a true critical point and will be 
suppressed to zero in the thermodynamic limit. 

Another finite-size effect is shown in Fig|Sl There are 
two obviously close apexes for both L = 10 and L — 12. 
For L = 8, there are two close apexes for smaller value of 
U. Actually, the first apex is not a critical point. It is also 
caused by the crossover between two PSs. Take L = 10 as 
an example, for V < —2.8, the ground state is dominated 
by configuration PS(d), while for -2.8 <V< -2.2, it is 
dominated by configuration PS (a). At first glance, this 
result seems unreasonable, because when V > —2.2, the 
ground state is dominated by SDW. Why the ground 
state tends to be doubly occupied before it becomes a 
SDW phase? A second order perturbation analysis will 
help us to understand the complication here. In the 
strong couphng limit, \t/U\ <C 1 and/or \t/V\ <C 1, we 
have the energy for the configuration PS(d) and PS(a), 
respectively, 5U+16V+4t^/V and 'iU+16V+{3/2)t^/V. 
For small \U\, the system prefers to be phase separated 
due to attractive Coulomb potential V, and the hopping 
process prefers the configuration PS (a). Thus the effect 
of hopping process is stronger than that of the on-site 
Coulomb interaction when the absolute values of V de- 
creases, and the dominating configuration changes from 
PS(d) to PS(a). Similar to the argument for Fig. |31 
this unexpected apex is also a finite size effect. Because 
when the block size is fixed, we vary system size to oo, the 
change from one pair of separate t and J, to a doubly oc- 
cupied will cause no weight change in the reduced-density 
matrix. 



FIG. 6: The block-block entanglement change with different 
values of U for fixed block size I = 4 and various system size. 
Here (7 = 2. 

VI. DISCUSSIONS 

An important observation in our study is that in order 
to reveal phase transition by the block-block entangle- 
ment, the choice of block type is essential. Sometimes, 
the block-block entanglement could be a smoothly con- 
tinuous function and shows no obvious structure at the 
transition point. As a result, finding the phase bound- 
aries of some transitions is not an easy job, when the 
block-block entanglement at the critical point behaves 
neither singular nor extremum. Extremal values or dis- 
continuities shown at some parameters is only suggestive. 
One must do more detailed analysis to identity where 
there exists true QPT and its nature. A typical example 
is the transition from SS to CDW, as shown in Fig. |21 
When such situation arises, one may ask oneself these 
questions: (i) What is the behavior of the block-block 
entanglement in the critical region, i.e. its derivatives? 
(ii) How about choose different block, such as the "block" 
in the momentum space instead of in real space due to 
off-diagonal-long-range order? 

Regarding the first question, Osterloh et.al}^ found 
out that in the one-dimensional transverse-field Ising 
model, the first derivative of the concurrence with re- 
spect to the coupling diverges at the critical point, even 
though the concurrence itself does not show maximum 
value at that point. This approach was also applied to 
other quantum spin modelsi^ and they found similar re- 
sults. In Fig. [7| we show the first derivative of the block- 
block entanglement with respect to the nearest-neighbor 
interaction V: dE/dV. At the point V = -0.2, the 
derivative of the block-block entanglement achieves a ex- 
treme value, which is very close to the phase transition 
point for the infinite chain a.tV = i). Therefore, besides 
its extremum behavior in the QPT process, the deriva- 
tive of the entanglement may show extremum at critical 
point. 

Regarding the second question, we choose the "block" 
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FIG. 8: The block-block entanglement in momentum space 
as a function of V . Here L = 8, Z = 4, and U = —2 

in the momentum space. The choice came from the na- 
ture of superconducting paring, for the order parameter 
is 

i 

thus, in the momentum space, one has, for singlet super- 
conducting phase, 

Ax + = — ^cosfcx4_-fC^^fc_^ (9) 
fe 

while for triplet superconducting phase, 

- ^^sin/cx4_^cLfe^^. (10) 

k 

So if we choose the block in the momentum space, i.e. 
±fc, the creating operator of cooper pairs, which is indeed 
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a superposition of the superconducting order parameter, 
is naturally included in the reduced density matrix. 

We compute the block-block entanglement in the mo- 
mentum space for the L — 8 system. Using the anti- 
periodic boundary conditions, the bases we choose are 
K = ±l/87r,±3/87r,±5/87r,±7/87r. The results of the 
block-block entanglement between K = ±l/87r and the 
rest of the momenta are shown in Fig. |S1 From the fig- 
ure we observe that the block-block entanglement in the 
momentum space shows a minimum value at the point 
V = —0.2. This minimum value may suggest a change 
of phase from SS to CDW. Surprisingly, this result is 
consistent with the result of the first derivative of the 
block-block entanglement in real space. However, the 
entanglement in the momentum space cannot be used to 
identify phase separation in the cross section. This is 
understandable, because phase separation occurs in real 
space. 

VII. SUMMARY AND ACKNOWLEDGEMENT 

In summary, we have studied the block-block entan- 
glement in the ground state of half-filled one dimensional 
extended Hubbard model. We found that it can iden- 
tify the main phases of EHM on the U — V plane, SDW, 
CDW, PS and even SS. Obviously, since the block-block 
entanglement with block size larger than one comprises 
nonlocal correlation, it provides more information than 
the local entanglement. Then a richer structure has been 
obtained from its contour map on the U ~V plane. More- 
over, the scaling analysis based on the block size implied 
that the log2(/) behavior of the block-block entanglement 
may also exist in the noncritical region, such as PS. This 



result is quite different from the previous studies of spin 
systems. We interpret it due to the high density-of-state 
closing to the ground state in this region. The number of 
different configurations taken part in the ground state is 
proportional to the system length. While in other non- 
critical regions, such as CDW, the ground state is nearly 
doubly degenerate and comprises two configurations, this 
fact restrict the increase of the block-block entanglement 
when the block size increases. By studying the depen- 
dence of the entanglement on system length, we also clar- 
ified its unusual behavior in the crossover from one PS 
to another and attribute it to the finite size effect. On 
the other hand, the derivatives of the block-block entan- 
glement and the block-block entanglement in momentum 
space also gave us some interesting results. 

However, it seems that the block-block entanglement 
still can not witness all phases of the system, such as 
the BOW and the TS. There are two possible reasons. 
One is that the system size we considered is not large 
enough, and the other is the limitation of the block- 
block entanglement itself. For the former, it is use- 
ful to apply other numerical method, such as density- 
matrix-renormalization-group, to this problem in the fu- 
ture studies. While for the later, it is necessary to search 
other kinds of proper way to quantify the entanglement in 
order to have a comprehensive understanding on the crit- 
ical behavior in the ground state of a many-body systems, 
just as we did in this work. Possible choices include sub- 
lattice entanglement suggested by Chen et al^^ and the 
fermion concurrence proposed by Deng and Gu"^*^. Obvi- 
ously, each choice has its own merits and limitations. 

This work is supported by the Earmarked Grant for 
Research from the Research Grants Council of HKSAR, 
China (Project CUHK 401504). 
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